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Transient dynamics in the Anderson-Holstein model with interfacial screening 
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We study the combined effects of electron-phonon coupling and dot-lead repulsion in the transport 
properties of the Anderson-Holstein model. We employ a recently proposed nonperturbative method 
to calculate the transient response of the system. By varying the initial conditions for the time 
propagation, we are able to disentangle two different dynamical processes, namely the local charge 
rearrangement due to the dot-lead contacting and the establishment of the noneqnilbrium many- 
body state due to the application of the external bias. According to the distinct initial contacting, 
the current can exhibit transient oscillations of different nature. These origin from tunneling events 
that involve virtual Franck-Condon excitations, or virtual transitions between the resonant level 
and the Fermi energy of the leads. 


PACS numbers: 71.38.-k,73.63.Kv, 81.07.Nb 

I. INTRODUCTION 

The interplay of vibrational and electronic de¬ 
grees of freedom during the charge tunneling through 
molecular junctions gives rise to striking nonequi¬ 
librium effects^. These include phonon-assisted 
tunneling^’^, hysteresis-induced bistability^^®, local 
heating^, molecular switching®’®, and negative differen¬ 
tial conductance®’^®’^^. Furthermore, if the electronic 
coupling to the vibrations is particularly strong a col¬ 
lective phenomenon known as Franck-Condon blockade 
(FCB) occurs®^. This manifests itself in a dramatic sup¬ 
pression of the tunneling current at low bias, accom¬ 
panied by the appearance of several vibrational side¬ 
bands inside the Coulomb-blockade diamonds, recently 
observed in carbon nanotubes quantum dots (QD)®®. 

From the theoretical side, the quantitative descrip¬ 
tion and understanding of the above phenomena repre¬ 
sents a challenging nonequilibrium problem. Many ob¬ 
served properties can be addressed within the Anderson- 
Holstein model®'^, that describes a single electronic level 
coupled linearly to a vibrational mode and to metallic 
electrodes. Despite its apparent simplicity, however, this 
model does not allow for an analytic solution and approx¬ 
imate approaches must be resorted. Only very recently, 
numerically exact methods to calculate its nonequilib¬ 
rium properties have been put forward®®^®®. Beside pro¬ 
viding a valuable validation of previous findings, the ex¬ 
act data have also revealed novel features, especially in 
the time domain, like the extraordinarily long-transient 
dynamics needed to reach the stationary state in the FCB 
regime®®. 

The inclusion of electron correlations widens the 
range of accessible phenomena, but at the same time 
complicates the theoretical treatment even further. 
While a considerable amount of work has been de¬ 
voted to study intra-molecule electron-electron (e-e) 
interactions^® (mainly to address the joint effects of 
vibrations and Coulomb blockade or Kondo-like corre¬ 


lations), the role of interfacial repulsion between the 
molecule and the leads has been explored only very 
recently®®. In that paper we developed a novel ap¬ 
proach to study the complex interplay between the lo¬ 
cal electron-phonon (e-p) coupling together with the 
molecule-lead e-e repulsion. We showed that the exci- 
ton formation at the molecule-lead interface improves 
significantly the polaron mobility thus competing with 
the FCB. As a consequence the FCB regime is dynami¬ 
cally established after a long-lasting sequence of blocking¬ 
deblocking events®®, characterized by exciton-enhanced 
current spikes. 

In this paper we extend our previous study by com¬ 
puting the time-evolution of QD density under the ap¬ 
plication of an external bias, and by investigating how 
the transient current is modified by changing the initial 
conditions for the time propagation. When possible, we 
also compare our results with the exact data, finding very 
good agreement. We show that different ways of connect¬ 
ing in time the QD and the leads can produce qualita¬ 
tive differences in the transient current, accompanied by 
a change of the nature of the dominant oscillations. 


II. MODEL AND FORMALISM 

The model we consider consists in a single-level QD 
attached to two semi-infinite one-dimensional noninter¬ 
acting wires. An electron occupying the level is coupled 
to a single vibrational mode located in the QD, and at the 
same time interacts with the electrons in the leads which 
are in the proximity of the QD. The spinless Hamiltonian 
describing this system is given by (in standard notation) 

OO 

H = -tw X + h.c.) -b Ti y^(dlid + h.c.) 

Oi^j — l OL 
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where a = L,R labels the Left and Right lead, and fid = 
dfd, fiao = dl^QdaO are the densities on the QD and on 
the first sites of the lead a respectively. The systeni is 
perturbed by an external bias given by Hy = 
with Na = Xx ^OLX the total number of particles in wire 
a and V = Vl — Vr the total applied voltage. 

In the following we consider the continuum version He 
of the above model since it allows for a nonperturba- 
tive treatment of the e-e interaction via the bosoniza- 
tion technique^^. To this end we first assume half-filled 
wide band leads with linear dispersion = vpk (with 
vp = 2tu,a the Fermi velocity and a the lattice spacing) 
and constant tunneling amplitude F = 27rT;^Xfc^(^ ~ 
efc) = 2T^/tw, and then we unfold the left and right 
leads^^’^^. In this way the first term of Eq. (I) takes 
the Dirac-like form — Xa / dxipl^{x)dx'il’a{x), where 
the electron field operator ipa{x) describes an electron 
at position x in the (chiral) lead a moving with velocity 
vp. The rest of the continuum model is simply obtained 
by replacing dax —t '4’a(x), Xx f dx, and by rescal¬ 
ing the model parameters according to Ti ^ ti = 2^/aTi 
and U ^ u = Aall. We then bosonize the electron op¬ 
erators as^^ ipa{x) = with boson field 

^a{x) = iaJ2q>oCg(biqe~'^°“^^ - h.c.) - ^xNajC, and 
ria the anticommuting Klein factor. In the mode expan- 

_av^ 

sion of the boson field it holds (q = ^ ^ 2 tq ’ dl the 
length of the system. The electron density in the leads 
takes the form ha(x) = —dx4>a(x)/y/p^ and hence, up to 
an irrelevant term^^, the bosonized continuum Hamilto¬ 
nian reads^^ 


evaluated in a; = 0 , with the effective interactions Wrr = 
WpL = 1 - u/{2pvp) and Wrl = Wrr = -u/{2pvp), 

\ 2 r> — CIQ 

and renormalized energy level Ed = ^d— ^ — u Xq ^vpC ■ 
In the new basis we have a noninteracting QD coupled 
to effective leads (bosonic baths) whose ground state for 
ti = 0 is l^'o) = |0p)(8)na<j lOa?)’ where |0p) and |0ag) are 
the vacua of the boson operators d and bag respectively; 
the tunneling coupling occurs via the correlated-polaron 
operator /. 


III. EQUATION OF MOTION 

The great advantage of casting the Hamiltonian in 
the form of Eq. (3) is the possibility of writing an 
approximate equation of motion for the nonequilibrium 
QD Green’s function, solely in presence of a correlated- 
polaron embedding self-energy. The latter accounts for 
the presence of the biased leads as well as for the e-p 
and e-p interactions in a nonperturbative way. We de¬ 
fine the QD Green’s function on the Keldysh contour as 
G{z,z') = j{Td{z)d\z')), where T is the contour order¬ 
ing, operators are in the Heisenberg picture with respect 
to H'e-f Hr (the bias perturbation does not change after 
the transformation); the average is taken over the un¬ 
contacted ground state |d>o) G \nd), \nd) being the state 
of the QD with density Ud, with Ud ranging from 0 to 
1. In Ref. 21 we proposed a controlled approximation 
scheme in order to derive a closed (and numerically solv¬ 
able) equation of motion for G, that reads 


He 


^Pdbiqb aq + Edfid + uJod)d 

a,q'>0 

+ h.c. 




nd 


A(d’^ -|- d) — M 


CqQ 


a,q>0 


Kq + ^aq) 
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Next we perform a multi-boson Lang-Firsov trans¬ 
formation to (formally) eliminate the e-p and e- 
p coupling appearing in the last line of the above 
equation^^’^"^. The unitary operator U = exp[—^(d’^' — 

d) -I- 2.y/TTuY,aq Sf^^aq ~ baq)]nd transforms the contin¬ 
uum Hamiltonian into id' = W Held with (from now on 
all sums are over g > 0 ) 


H'e = '^vpqb'lqbaq-fuJoa^d-fednd + ti'Y^ ^flgd-f h.c. . 

aq a 

( 3 ) 

In the transformed Hamiltonian it appears the renormal¬ 
ized fermion field 

( 4 ) 


{idx - £d)G{z, z') - dz'^ Ea(z, z)G{z, z') = S(z, z'), 

dl a 

( 5 ) 

where Sq,(z, z') is the correlated-polaron embedding self¬ 
energy. The real-time Keldysh components of S can be 
evaluated exactly using again the bosonization method^^ 
and read 


s|(t,o = ± 


ivpTe ^ i\q,^(t')-(pa(t)\ 


age 


±-iu>Q (t — t'} 


[a =F ivp{t — t')\d 


,( 6 ) 


with adimensional e-p coupling g = (A/wq)^, interac¬ 
tion dependent exponent jd = 1 -\- and phase 

(pa{t) = J*diVa{t). The power-law reflects the collective 
excitonic response of the lead electrons to the attractive 
potential due to the creation of a hole in the QD^®. In 
the noncorrelated case we have fd = 1, while electron cor¬ 
relations produce /3 < 1 ; the smaller the exponent is, the 
stronger is the exciton effect. 

The integral in Eq. (5) runs over the Keldysh con¬ 
tour 7 , and using the Langreth rules^® it is converted 
into a coupled system of Kadanoff-Baym equations®®’®^ 
(KBE) which we solve numerically. Once the Keldysh 
components of G(z, z') are known, the time-dependent 
QD density is calculated as n{t) = —iG’^{t,t), while the 
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FIG. 1: TD density in the partitioned approach for [7 = 0 
with initial QD occupancy = 0 (red) and rid = 1 (green). 
Exact data from Ref. 16 are also displayed (circles). The 
rest of parameters are A = 16, ajo = 8, E = 26, Id = —10, 
vf/cl = 100. Units: F for energies and F”^ for times. The 
insets displays the TD result for a longer propagation time, 
in order to appreciate that the two densities reach the same 
steady-state value. 

transient current flowing through the QD and the a lead 
can be evaluated according to 

Ia{z) = / dzTia{z,z)G{z, z) + h..c. (7) 

Jy 

IV. INITIAL CONDITIONS 

Before presenting the explicit numerical results, a brief 
discussion about the initial conditions is in order. To 
solve practically Eq. (5), one has to set the initial value 
G'<(0,0) = —iud, that corresponds to (apart from the 
—i factor) the initial density in the QD. This means that 
our propagation scheme assumes the system initially un¬ 
contacted {ti = 0, i.e. the leads and the QD are in their 
ground states), and the contacts and bias are switched 
at time t = 0. Therefore the subsequent transient regime 
accounts for two different dynamical processes: (i) the 
charge rearrangement due to the QD-lead contacting 
(creation of Friedel-like oscillations in the leads, etc), and 
(ii) the establishment of a genuine nonequilbrium many- 
body state due to external bias (the rise of the current 
towards its steady-state value, etc). This scheme corre¬ 
sponds to the so-called partitioned approach. However, in 
a more realistic situation, the bias is switched when the 
system is already contacted and in equilibrium (i.e. the 
leads and the QD are in the ground state with ti yf 0). 
This is the so-called partition-free approach^^’^^. Here 
the transient dynamics can be very different from the 
partitioned case, since the two processes (i) and (ii) de¬ 
scribed above are not superimposed^"'’’^®’^®’^®’^®’^®. In our 
scheme we can numerically simulate the partition-free ap- 
praoch, since the bias function Va{t) (appearing in the 
phase ifa in Eq. (6)) is completely arbitrary. In practice 
we consider a step-like bias function Va{t) = Va9{tth), 
that corresponds to take the system initially uncontacted 



FIG. 2: TD relaxation towards the equilibrium density for 
U = 0 (black) and U = (red), with initial QD occupancy 
Ud — 0.9 and bias V = 0. The rest of parameters and units 
are the same as in Fig. 1. 


with a given at time t = 0, let the system thermalize 
(dynamics (i)) till a time tth at which no current flows 
across the links, and then we switch the bias perturbation 
on (dynamics (ii))"^®. 


V. TRANSIENT DENSITY 

In Fig. 1 we assess the accuracy of the proposed ap¬ 
proach by comparing our results against exact data avail- 
abale in the literature for {7 = 0, and obtained within the 
partitioned scheme^®. It appears that the agreement is 
exceptionally good for initial density = 1, while for 
nd = 0 we predict a slower raise of the density towards 
its steady-state value. We recall, however, that in this 
case our results improve the state-of-the-art^®. In the 
inset we show the TD density for a longer propagation 
time (not within reach of current numerical techniques), 
in order to appreciate that the two densities reach the 
same steady-state value, as it should be. 

We now study the effects induced by the QD-lead re¬ 
pulsion U. Fig. 2 displays the relaxation of the QD 
density from the uncontacted value (ud = 0.9) to the 
equilibrium value, after that the system has been con¬ 
tacted without bias. We see that the effect of the screen¬ 
ing interaction is twofold: enhance the asymptotic value 
of n(t), and speed up the relaxation time. The first effect 
origins from the fact that a finite electron density in the 
QD induces a charge depletion in the portion of the leads 
which are in the proximity of the interface; part of the 
repelled charge, in turn, migrates towards the QD thus 
enhancing its population. The second, instead, is a jam¬ 
ming effect reminiscent of the one found in Ref. 24, that 
is due to the dynamical screening of the QD charge and 
tends to stabilize faster the value of n{t). 
























4 



FIG. 3: TD current within the partitioned scheme for different 
U, initial QD occupancy rid = 1, and bias voltage V = SF. 
For U — 0, exact data from Ref.^® are also displayed (circles). 
The rest of parameters and units are the same as in Fig. 1. 


VI. TRANSIENT CURRENT 

In this section we study the TD current I{t) = + 

Iii{t)]/2 flowing under the application of the external 
bias. As in the case of the density, the approach is first 
validated by comparing our results with the exact re¬ 
sults recently obtained for [7 = 0 within diagrammatic 
Monte Carlo simulations^®. In Fig. 3 a remarkable agree¬ 
ment between the two approaches (within the partitioned 
scheme) can be appreciated. In particular our method ef¬ 
ficiently reproduces the very peculiar transient behavior 
of I{t) before the steady-state is reached. In the absence 
of e-e interactions the time-dependent current displays 
quasi-stationary plateaus where almost no electron tun¬ 
nels across the junction. At times tn = 2mr/u!o a de¬ 
blocking effect occurs, and the current exhibits narrow 
bumps, signaling a sudden electron flow. When the e- 
e interaction is considered, we observe a significant en¬ 
hancement of the transport properties, characterized by 
larger current spikes. As we show below, the physical in¬ 
terpretation of such striking transient dynamics provides 
a clue to understand how the FCB regime is dynamically 
established, and how e-e interaction modify the FCB sce¬ 
nario. At t = tn an electron occupying the QD is in 
the polaron ground-state with phonon cloud centered at 
X ^ nX (with n ~ 1). At this time the polaron tunnels 
to the lead, causing a displacement of the oscillator to 
cc —>■ 0 (since now n = 0). At this point the polaron can¬ 
not hop back to the QD since the overlap between the 
two shifted oscillator wavefunctions is negligible. Only 
after a vibrational period 2'k/ujo the overlap returns to 
be large, and the polaron can hop back to the QD. Let 
us now consider the effects of e-e interaction. The mech¬ 
anism described above is modified as follows: If at time 
t = tn ane electron is on the QD, the electron density 
diminishes at the site of the lead boundary, thus over¬ 
coming the hopping suppression due to the Pauli prin¬ 
ciple and enhancing the effective tunneling rate^^’"^^’^^. 
Similarly, at time t = tn+i the electron can easily tunnel 
back to the QD, being attracted by the hole previously 




FIG. 4: Comparison between TD currents within the partion- 
free scheme (black) vs partitioned scheme with = 0 (blue) 
and rid = 1 (red). [7 = 0 in the upper panel and U — t^ in the 
lower panel. The rest of parameters and units are the same 
as in Fig. 1. The insets display the TD result for a longer 
propagation time. The TD curve of the partition-free case 
has been shifted of tth = 100. 


left. This explains in a transparent way the [7-induced 
enhancement of the current spikes observed in Fig 3. 

We now focus on the effects of the different initial 
conditions. To this end we compare the results ob¬ 
tained within the partitioned scheme vs the ones ob¬ 
tained within partition-free scheme. In Fig. 4 we plot 
the TD current in the two schemes without (upper panel) 
and with (lower panel) screening interaction [7. For the 
partitioned scheme we employ two different initial val¬ 
ues of the QD density, namely = 1 (red curve) and 
rid = 0 (blue curve). We recall that in order to display 
the partition-free curves together with the partitioned 
ones, we shift the former of tth- We observe that the 
three currents correctly reach the same steady-state and 
do it via similar long-lasting sequences of current spikes, 
(typical of the FCB regime), and, as expected, the cur¬ 
rent in presence of U is enhanced with respect to that 
calculated at [7 = 0 (Coulomb deblocking). In this case 
the partitioning effects reflect in a quantitative change of 
the early transient current (see the black curves vs the 
red and blue ones), even though the qualitative FCB-like 
transient behavior does not change. 

In order to magnify the partitioning effects, we focus 
on the perfect resonant case, obtained by setting id = 0. 
Here (see Fig. 5) the partition-free transient current 
displays qualitative differences. If we use Ud = 1/2 for 
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FIG. 5: Comparison between TD currents in the resonant con¬ 
dition Sd = 0 within the partion-free scheme (dashed black) vs 
partitioned scheme with = 0.5 (blue). 17 = 0 in the upper 
panel and U = tm in the lower panel. The rest of parameters 
and units are the same as in Fig. 1. The insets display the 
TD result for a longer propagation time. The TD curve of 
the partition-free case has been shifted of tth = 100. 


TD density is identical in the two schemes, the TD cur¬ 
rents are remarkably different. In particular the FCB 
pattern tends to disappear in the partition-free case, and 
the transient current has a smooth oscillating behavior 
with dominant frequency F/2, i.e the energy difference 
between the molecular level and the Fermi energy of the 
leads. Here the interaction U (lower panel of Fig. 5) am¬ 
plifies these oscillations, thus enhancing further the dif¬ 
ference between the two currents. This behavior is con¬ 
firmed by decreasing the bias from V = 26r to F = SF, 
as shown in Fig. 6. It is clear that the transient fre¬ 
quency of the partition-free current reduces according to 
the smaller value of V, whereas the spikes of the parti¬ 
tioned current continue to appear with periodicity given 
by the phonon frequency uq. 

This qualitative change can be interpreted as follows: 
At the resonance (id = 0) and if the leads are initially 
relaxed around the QD, the dressed tunneling (that in¬ 
volves the excitation of a phonons with energy uiq) has 
the same probability as the bare tunneling (that involves 
virtual transition between the QD and the Fermi level 
of the leads). In the partitioned case, instead, the bare 
tunneling is suppressed because the electron charge in 
the leads is optimally rearranged in the proximity of the 
QD. 


VII. CONCLUSIONS 



FIG. 6: Comparison between TD currents in the resonant con¬ 
dition £d = 0 within the partion-free scheme (dashed black) 
vs partitioned scheme with n^j = 0.5 (blue). U = Q and bias 
V = 5r. The rest of parameters and units are the same as 
in Fig. 1. The TD curve of the partition-free case has been 
shifted of tth = 100. 

partitioned scheme, the TD density is pinned in both 
schemes at the constant value n(t) = 1/2 at each time 
(not shown), due to symmetry constraints. Despite the 


We presented a systematic study of the time-dependent 
transport properties of the Anderson-Holstein model in 
presence of dot-lead screening interaction. Thanks to 
the analytic expression of the correlated-polaron embed¬ 
ding self-energy we gained a clear understanding of how 
the two interactions combine together. The validity of 
the approach was also corroborated by comparing our re¬ 
sults against exact data available in the literature. Our 
approximated scheme allows for a calculation of the time- 
dependent density that improves the current state-of-the- 
art, and at the same incorporates the screening effects 
in a physically correct way. The transient behavior of 
the current was carefully analyzed by considering differ¬ 
ent initial contacting of the leads. We showed that at 
early times the current can exhibit qualitative different 
behaviors, depending if the electron liquid in the leads is 
relaxed or not around the QD before the switching of the 
bias. At resonance, we found that the partition-free cur¬ 
rent displays coherent oscillations with frequency equal 
to the applied bias, whereas the partitioned current has 
a periodicity dictated by the phonon frequency. 

We acknowledge funding by MIUR FIRB grant No. 
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